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We derive the leading non-global logarithms (NGLs) of ratios of jet masses mi, 2 and a jet energy 
veto A due to soft gluons splitting into regions in and out of jets. Such NGLs appear in any 
exclusive jet cross section with multiple jet measurements or with a veto imposed on additional 
jets. Here, we consider back-to-back jets of radius R produced in e+e~ collisions, found with a cone 
or recombination algorithm. The leading NGLs are of the form ln^(A/mi,2) or In'^C^i-i/'^z). 
Their coefficients depend both on the algorithm and on R. We consider cone, kx, anti-kx, and 
Cambridge-Aachen algorithms. In addition to determining the full algorithmic and R dependence 
of the leading NGLs, we derive new relations among their coefficients. We also derive to all orders 
in Ua a factorized form for the soft function S{kL, ku, A) in the cross section a{mi,m2, A) in which 
dependence on each of the global logs of n/kL, fi/kn and /x/A determined by the renormalization 
group are separated from one another and from the non-global logs. The same kind of soft function, 
its associated non-global structure, and the algorithmic dependence we derive here will also arise in 
exclusive jet cross sections at hadron colliders, and must be understood and brought under control 
to achieve precise theoretical predictions. 



I. INTRODUCTION 

Observables of varying exclusivity can be used to probe 
the jetlike structure of final states in high energy colli- 
sions [2 12] or even the substructure of the jets them- 
selves [3]. While more exclusivity reveals more infor- 
mation about structure, it also introduces dependence 
on additional scales, ratios of which induces potentially 
large logarithms in perturbative expansions. 

Non-global observables [H [S] are those for which soft 
radiation in sharply divided regions of phase space are 
probed with different measures. For example, measur- 
ing separate masses toi,2 of two back-to-back jets pro- 
duced in e+e~ collisions generates non-global logarithms 
(NGLs) of nil/ 1712- Measuring the total invariant mass 
of a two-jet-like final state, = -I- rn^, however, 
generates no NGLs since the soft radiation everywhere is 
probed "globally" with a single scale m. Global logs like 
those of m/Q in such a cross section can be resummed 
using well-known methods (e.g. [6]). The most power- 
ful of these are based on the renormalization group (RG) 
evolution of hard, jet and soft functions in a factorization 
theorem for the global observable in perturbative QCD 
[li ,8. or soft-coUinear effective theory (SCET) [9Hl3]. 

To be more precise, for the example of the hemisphere 
dijet mass distribution, a factorization theorem for the 
distribution in 2 takes the form [TH [15] 

xJn{ml-Qkii,fi)S{kL,kB,,fJ.)-\ , 

where k^ fi measure the small light-cone components 
n-k,n-k of the total momenta in each hemisphere, where 
n,n = (1, ±z) with +z along the jet axis in the right (R) 



hemisphere. The soft function S can be expressed as a 
convolution of perturbative and nonperturbative pieces 
|16) . but we only consider its perturbative component 
here. The hard function H depends only on logs of fi/Q 
and the jet functions J1.2 depend on logs of /i/TOi,2, but 
the soft function depends on logs of ^/k^, fJ-fk^, and 
kh/kR. The /^-dependent "global" logs can be resummed 
by RG evolution, but the logs of kL/ka are non-global 
and not resummed by the ordinary RG. In other words, 
logs of the ratio of any single soft scale to the hard scale 
or either jet scale can be resummed by running between 
those scales. But logs of ratios of soft scales among them- 
selves cannot be resummed by using a framework that 
only contains one soft mode and thus one soft scale to or 
from which RG evolution can be performed. 

Another case in which multiple soft scales appear is 
when two jets i, j in a multijet event come close together. 
This introduces the soft scales rri^/mij and mj/mij, as- 
sociated with the "fat" dijet, in to the cross section, in 
addition to the usual soft scales, mf/Q, associated with 
individual jets. Logs induced by these additional soft 
scales can be summed using the effective theory SCET+, 
an extension of SCETi that contains a "collinear-soft" 
(csoft) mode [T7]. This was the first factorization of 
a multi-scale soft function that allows for complete re- 
summation of all the resulting large logarithms in the 
cross section. These logs should be distinguished from 
non-global logs, which come from making measurements 
in different regions of phase space and are ratios of soft 
scales (e.g., nii/mj). 

NGLs were first recognized in [1] in e^e^ dijet event 
shape distributions in which the mass of only one hemi- 
sphere jet PR — mif /Q is measured while being inclusive 
in the other hemisphere, and in a larger class of event 
shapes in deep inelastic scattering in . Subsequently, 
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Refs. [ini 120] studied NGLs of K/Q in cross sections 
vetoing radiation with total energy greater than A in an- 
gular regions outside of found jets. Though a hard scale 
Q appears in these ratios, we found in [2r that the NGLs 
still arise from considering both scales in the ratio to be 
soft and later taking one of them to Q in an inclusive 
limit. 

In j21| we made progress in understanding the ori- 
gin of NGLs in effective field theory. We considered the 
factorized dijet invariant mass distribution (7(7711,7712) in 
e+e^ collisions producing back-to-back jets, and calcu- 
lated to ©(a^), as also in [22], the hemisphere soft func- 
tion S{kL,kji). These calculations clarified the origin of 
NGLs in an EFT framework as the dependence of a soft 
function on ratios of multiple soft scales, and revealed 
new subleading (single) NGLs and non-logarithmic non- 
global functions. 

These NGLs are organized into a multiplicative factor 
entering the total cross section, with the leading NGLs 
taking the generic form 

2 

5NG(/ii/A^2) = 1 - -^CfCaS2 In^ ^ + • • • . (2) 

Here /ii_2 are the scales at which soft radiation is probed 
in different sharply-divided regions. For the hemisphere 
mass distribution /ii^2 = "m-i 2/Q ^-nd S'2 = 7r^/3. For 
the PR distribution, fii = Qpr while P2 = Q due to 
total inclusivity in one hemisphere. The coefficient S'2 
is a geometric measure of the region into which the two 
soft gluons contributing to a NGL can go. The fact that 
it varies with the size of this region is due to the NGL 
arising from a purely soft divergence of QCD. Techniques 
to resum NGLs using numerical fits in the large- A^c limit 
of QCD were introduced by |1] , but analytic resummation 
of NGLs in real-world QCD remains an open problem. 

In this work we seek to extend the intuition gained in 
PT] by studying a more exclusive set of cross sections. 
We study non-global properties of an exclusive jet cross 
section (j{mi, m2, A), where the invariant masses 7771 and 
7M2 of two jets of size R produced in an e+e^ collision 
at center-of-mass energy Q are measured, with a veto A 
on the energy of additional jets. We consider finding the 
jets using various algorithms — cone, anti-kx, Cambridge- 
Aachen, and kT [23ti28] . We will find that NGLs of 
the ratio of the jet veto and the jet masses K/mi^2 
are present, in addition to NGLs of the ratio of masses 
mi/m2- We calculate the coefficients only of leading dou- 
ble NGLs In^ (1^1/1^2) in this paper. The relevant scales 
for this observable are shown in Fig.[l]for a particular hi- 
erarchy of 7771^2 and A, however our results are valid for 
any choice such that Q ^ 7^71.2 ^ mf 2/Q1 A. 

In |21J, we discovered that at 0{a1) NGLs of two soft 
scales /7i.2 can be constructed from separate pieces de- 
pendent on the ratio of the factorization scale /i to one 
physical scale at a time. Namely, the region of phase 
space where one of the soft gluons enters the region sen- 
sitive to the scale pi and the other enters the region 
sensitive to p2 generates the double log In^ /i^/ (/ii/i2). 




Right jet scale - "^2 



f /i^ = mj/Q 



Soft scales = A 



//f = ml/Q 



FIG. 1: The relevant scales in the exclusive jet mass cross 
section with an energy veto, A outside of the jets is shown 
for a particular choice of the hierarchy -012 <C AQ <C ml that 
gives rise to large non-global logs. Our results apply to any 
choice of 7^,1,2 and A that satisfies Q ^ mi, 2 3> m^ 2/Q, ^, 
which maintains the separation between hard, jet and soft 
scales. 



while the regions where soft gluons enter only region 1 or 
only region 2 generate \n^{fi/fii) and (fi/ IJ.2) ■ In 

|21) we derived from RG invariance of the cross section 
and IR safety of the soft function that the coefficients 
of these logs are constrained so that the /i-dependence 
cancels, but an NGL ln^(/ii//i2) is left over. Analo- 
gously for cr(777i, 7772, A), the three soft phase space re- 
gions that give rise to the NGLs at ©(a^) are shown 
in Fig. [2] Each configuration contributes logarithms of 
fi over a single scale, the "in-out" regions contributing 
logs In^ /Lt^/ (A 7771^2), and the "in-in" region contribut- 
ing logs In^ /x^ 7(77717772). These combine with single- 
region contributions to give NGLs of A/n7i_2 with coeffi- 
cients /oL,OR and of mi/m2 with coefficient /lr. These 
coefficients give the geometric factor S'2 in Eq. ([2]). IR 
safety and RG invariance will allow us to derive addi- 
tional strong relations among these different coefficients. 
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The division of the NGL into separately /i-dependent 
pieces opens up the possibihty to use the renormahzation 
group to sum NGLs, ahhough this has yet to be carried 
out expHcitly. 

It is worth noting that our calculation of the coeffi- 
cients /oL,OR,LR of the leading NGLs for a two-jet con- 
figuration applies to other measurements of the soft ra- 
diation in and out of the jets as well. Different choices of 
observable change the arguments of the NGLs, but their 
coefficients are related to the geometry of the configura- 
tion given in Fig.[2]and will be given by one of /ol,or,lr- 

While the presence of NGLs has previously been rec- 
ognized in jet mass observables similar to the one we 
consider here, we derive the full algorithmic and R de- 
pendence of the leading NGLs involving both jet masses 
and vetoes here for the first time. Refs. [SI [12] derived 
the NGL of A/Q in cross sections with a veto region 
outside jets defined with fixed cones (but with masses 
unmeasured). Ref. [2^ studied (and [29] calculated more 
precisely) the effect of clustering soft gluons with a fixed 
R = 1.0 with the kx algorithm, with a variable rapidity 
gap Ar] describing the size of the region in which radia- 
tion is vetoed. Ref. [30] considered a cross section with 
two cone jets of radius R <^ 1, with only one jet's mass 
m being measured and a veto A on radiation outside the 
jets, and calculated the NGLs of to/ A due to the mea- 
sured jet and A/Q due to the inclusive jet. Ref. was 
the first to study NGLs in an observable probing jet sub- 
structure, the mass drop in filtered subjets produced by 
decays of boosted Higgs bosons [32]. Our work consol- 
idates and extends many of these results by calculating 
NGLs of A/mi^2 and toi/to2 with two measured jets for 
four different cone and recombination algorithms, as a 
function of arbitrary jet size 0<i?<7r/2^, and deriv- 
ing new relations among the coefficients of the different 
NGLs. 

For the kx algorithm, it was pointed out in [5H1|3S] and 
elaborated in [30] that the effects of soft gluon clustering 
also affect the independent emission contributions to the 
types of observables mentioned above. The factorization 
theorems and resummed predictions we consider below 
should be modified to include such effects for clustering 
algorithms. We leave this outside the scope of our present 
work, focusing just on how NGLs affect such predictions. 

We gain some intuition from our investigation that is 
not necessarily novel, but is often not appreciated, and 
that we hope is helpful to clarify. First, the presence of 
multiple soft scales is enough to induce NGLs, regardless 
of their ordering. Second, we emphasize that NGLs arise 
not only from soft gluons splitting right along jet bound- 
aries and entering just inside the respective separate re- 
gions, but from the entire angular region of phase space 
into which the soft gluons can enter. The numerically 



largest contribution comes from the two gluons being 
close to each other near the boundary, but the enhance- 
ment is not parametrically large. Thus NGLs cannot 
be avoided simply by erasing dependence on the bound- 
ary region in constructing a jet observable, although they 
can be somewhat reduced. Probing separate regions with 
separate soft scales is enough to induce NGLs [50]. 



While our results here are derived for e~^e~ collisions, 
the methods and lessons are directly applicable to exclu- 
sive jet cross sections measured at hadron colliders such 
as the LHC. For instance, distributions in multiple jet 
masses, jet shapes or event shapes such as A^-jettiness 
[36l l37] with different values of the measure on the N 
jets will contain NGLs of mi/nij or ti/t2. Exclusive jet 
cross sections defined with explicit vetoes on the pr of 
additional jets induce NGLs of p™* over the relevant hard 
scale. For example, vetoing jets in searches for Higgs to 
WW — >■ iviv with zero jets will induce NGLs oip^^ /tuh- 
It would thus be wise either to calculate and control these 
NGLs or to use methods to veto jets that avoid NGLs, 
such as beam thrust (0-jettiness) based vetoes [5Slf551f55] . 



A final consequence of our results regards the claim of 
the recent Ref. f40| by Kelley, Schwartz, and Zhu (KSZ) 
that the cross section cr(/9. A), where p — {m\ +m'^)/Q'^, 
contains no logs of A/{Qp), and that therefore a factor- 
ization of the form a{p,A) = crin(p)crout(A) holds to all 
orders in a^, at least in the regime A < Qp <^ QR <^ 1^. 
This conclusion is not consistent with our calculations, 
which show that NGLs of A/{Qp) are in fact present for 
any value of this ratio. We verified our prediction for the 
coefficient of the NGL by comparing to the predictions 
of the numerical Monte Carlo EVENT2 [ITIH^. 



In Sec. |TT] we review the factorization theorem for 
the cross section cr(TOi,TO2,A) and derive to all orders 
in as a correctly factorized form for the soft function 



S{kL,kji, A) that appears therein. In Sec. Ill we derive 



new generic relations among coefficient s of the different 
NGLs appearing in S{ki^, A;^, A). In Sec. IV we derive the 
leading NGLs appearing in S{kL , kfj, A) , including the 
full algorithmic and R dependence. In Sec. |V] we com- 
pare the predictions of cr{p, A) with and without the pre- 
dicted NGL of Qp/A to those of EVENT2 at 0{al), and 
confirm the presence and predicted sizes of the NGLs. 
In Sec. |VI| we conclude. In two Appendices, we provide 
the ingredients necessary to construct the global logs in 
ct(toi,TO2,A) from the RG, and the Feynman diagram 
amplitudes necessary to calculate the leading NGLs. 



^ We do require _R to be large enough for the observable we consider 
to factorize, namely mi, 2 Qtan(R/2) 1331 134| . 



^ We note that this claim appeared in version 1 of [40) . and has 
since been retracted in later versions. 



4 




In-Out Configuration, /or In-Out Configuration, /ql In-In Configuration, /lr 



[a) (&) (c) 

FIG. 2: The three basic soft gluon configurations that we use to calculate the leading NGLs. At 0{a^), the coefficient /or 
(/ol) receives contributions from one gluon in the right (left) jet and one gluon out of both jets, and the coefficient /lr receives 
contributions from one gluon in each jet. 



II. SOFT FUNCTION FOR TWO JET MASSES 
AND A VETO 

Factorization and resummation of (global) logarithms 
for exclusive jet cross sections defined with cone or re- 
combination algorithms with an energy veto outside the 
jets was first performed in [33l |34j. They imply that 
a {mi, 1712, A) factorizes in the form 

cT(mi,m2, A) = aoH{Q;^) J dkLdkRJn{ml - QkL] fi) 
X Jnimj - Qk^:, fJ,)S{kL, kii, A; R; fi) . (3) 

CTo is the Born cross section for e+e^ — > qq, Jn,n are jet 
functions, and S is the soft function. The scales that ap- 
pear in this factorization theorem are depicted in Fig. [T] 
We will also consider the cumulant distribution, defined 

by 

rQpi rQP2 /-A 

T,(pi,p2,A) = / dnii / dm2 / dA'a{mi,m2, A') 



which also factorizes in the form [21] 
E(pi,p2,A) = aoH{Q,fi) J dkLdk^JniQpi 



(4) 



(5) 



X JniQp2~kii, p)S''{kL, kn, A; R; p) , 



where 5^ is the cumulant soft function. The jet func- 
tions also depend on the jet size R and on the algorithm 
[3311311115, 44 but, as shown in these references, in the 
limit mi_2 <C (5tan(i?/2), the dependence on R is power 
suppressed, and Jn.n are the usual inclusive jet functions 
[^51 - H7| . We will work in this limit in what follows. 

The soft function S{kL, kn, A; R) arising in Eq. ([s]) was 
first defined and calculated to 0{as) in [33l [34] and is 
given by 



SikL,kn, A; i?) = ^ I T[Y^Y^] 

ISC ' 



X (5 A 



iex 



Xs 



(6) 



y,5(k 



iex 



ertn-kAsik 



iGX 



The theta functions 6in,out choose those soft particles i 
that end up inside one of the jets or outside both jets. 
Their precise form depends on the algorithm. In this 
definition, the energy outside the jets is fixed to be A, 
but we can integrate to obtain the cumulant which allows 
all energies up to A. 

Much of the structure of the soft function is determined 
by consistency of the factorization theorem in Eq. ([3| and 
the RG evolution of the hard, jet, and soft functions. We 
will argue its perturbative structure must take the form. 



S{kL,kR,A;fj.) = [SinikL; p)Sinik]i; p.)]Sout{A; p) 
(g) SjqcikL, fcfl;, A) , 



(7) 



where denotes a convolution of S'ng with S'in's in the 
variables k^n and with S'out in the variable A. The cu- 
mulant soft function S'^ behaves similarly. The pieces 
'S'in.out are determined by RG evolution and S'ng is not. 
'S'in,out depend individually on the scales fc^, ku, and A, 
while the S'ng has non-separable dependence on the ra- 
tios k^/ki^ and k^^j^/A. 

We can derive the form Eq. ([7| of S from RG invariance 
of the cross section cr, which is //-independent. Since 
the hard and jet functions (strictly speaking, its Laplace 
transform) have anomalous dimensions of the form 



fi—hiF = Tpln -I- jF , 



dp 



(8) 



where pp — Q for F — H and pp — Q{pi,2)^^^ for F = 
Ji^2, the soft function must have an anomalous dimension 



p^lnS^^Tsln^ 
dp 



In 1^ + 7s , 



(9) 



where Fs = —Th/2 — Tj and 7s = —Jh — 27,/. Notably, 
rS;7s are independent of R, and the scale A cannot ap- 
pear in Eq. ^ since the hard and jet functions know 
nothing about A. The pieces Th,j,s are proportional to 
the cusp anomalous dimension Fcusp to all orders in ag- 
in [33l [34], we calculated contributions to the soft 
anomalous dimension at 0{as) into pieces coming 
from a gluon inside a jet or outside the jets, finding 
d{liiS'')/d{lnp) = 7l + 7fi -f 7out, where 



1L,R 



Fein 



/ p tan ^ \ 
V kr.R ) 



(10) 
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with Ts the same as in Eq. (|9| and 7in = at 0{as 
Thus, 



7out = -2r5lntan' 



R 



2is - 27i. 



(11) 



which was verified by direct calculation to 0{as) in [331 
|34]. (Note 7s = also at 0{as).) 

The soft anomalous dimension can always be split up 



additively into pieces 7L.fl,7out of the form in Eqs. (10) 



and (111 to all orders in Ug. Finding operator or phase 



space definitions for soft functions which have these as 
their anomalous dimensions to all orders in as is another 
question. We found such definitions which work to 0{as) 
in ^33ti34j- Those definitions suggest that the pieces 7l,/j 
arise from the UV poles in soft function diagrams with 
all gluons inside the L,R jets, and 7out from the UV 
poles in diagrams with all gluons outside the jets. We 
will call the sum of all such diagrams Sl,r and ^out, 
respectively. Here, we only need these sums to be defined 
as those which have Eqs. ( 10 1 and ( 11 ) as their anomalous 
dimensions. It would take further work to show that 



their operator definitions in [33l |34| have precisely these 
anomalous dimensions to all orders, but these definitions 
are not necessary to derive the generic form Eq. ([7| of 
the perturbative soft function. 

The anomalous dimensions in Eqs. (10 1 and ( [TT| ) im- 
ply that the full soft function S'^{kL,k^^; R; fi) evolves 
under the RG according to 

S^ikL, kR, A; R; n) = S^ikL^kji, A; R\ ^o)Uout{R; Mi A*o) 
(8) [UL{kL\R\ Mi y'o)UR{kri,; R; ^, iiq)] , 

(12) 

where the evolution kernels Up — exp[fCi?(/i, /io) + 
LOpjlJ', Mo) 1^(m /M-f) ]7 with Kf,l^f defined in terms of 
rF,'yF in Eq. 



Splitting up the evolution kernels in Eq. (12 1 cleverly, 
we obtain 

S^{kL,kR,A;R;fi)=S''{kL,kR,A;R;no) (13) 
(g) C/out (i?; M, 2A) [Ul {kL]R; M, kL)UR{kR- R; fx, kR)] 
(g) Uout{R;2A, fio)[UL{kL; R; kL, fJ-o)UR{kR- R; kR, fio)] . 

The last two factors in brackets, finally, can be further 
split into 



UL{kL]R; kL,2K) ® UL{kL;R; 2A, ^lo) 

«) UR{kR; R- kR, 2A) ® UR{kR; R; 2A, /iq) • 



(14) 



Now all the /^Q-dependent factors in Eqs. (13) and (14) 



simply evolve S'^{kL, kR, A; R; fio) in Eq. ( 13 1 from /io to 



2A, so no Mo dependence actually remains. Then the soft 
function (and therefore S) automatically takes the 
factorized form of Eq. ([7| , with 

S[^ikL,^i) = UL{kL■,R■,^i,kL) (15a) 
SinikR, m) = UR{kR; R; m, fc^) (15b) 
5„^„t(A,/x) = C/out(i?;M,2A), (15c) 



and 

SI,G{kL,kR,K)=S'{kL,kR,K;R,2K) (16) 
X ULikL-, R; kL,2A)URikR; R; kR, 2A) . 

All logs in iS^ associated with the anomalous dimensions 
are now removed by evaluating S'' at fi = 2A in the first 
line and evolving between 2A and k^ R with the Ul.r 
kernels in the second line. The anomalous dimension of 
S in Eq. ^ then comes entirely from 5in,out , whose forms 
are now completely determined by the form of ^L,R,out 
and the RG. The remainder S'ng contains the dependence 
on the dimensionless ratios k^/kR and A/k^^R in S not 
constrained by the RG. 

It is notable that the form Eq. ([7| follows just from 
the form of the soft anomalous dimension Eq. (|9| that 
follows from RG invariance. It would take further work 
to give operator or phase space definitions of Si^ out that 
have precisely the anomalous dimensions Eqs. ( 10 1 and 



( 11 ) to all orders in as, but the generic form Eq. ^Tf does 
not depend on such a construction. 

Below we will give phase space definitions for S'iii.out,NG 
that work to at least ©(a^). We give results for the 0(0^) 
terms in i'in.out predicted by RG evolution in App.jBj and 
calculate the leading NGLs in S'ng in Sec. |IV[ 



III. 



STRUCTURE OF THE NON-GLOBAL 
DOUBLE LOGS 



While the renormalization group does not determine 
the NGLs in the soft function in Eq. ([3| , constraints from 
RG invariance can still be used to simplify the method 
to calculate them [21]. Here we use these constraints to 
derive relationships among the different "in-in" and "in- 
out" NGLs in the observables that we consider. 

To study the relationships between NGLs, we take the 
two jet radii to be different, Rl and Rr for the left and 
right jets respectively. While this can lead to ambigui- 
ties about which R to use in recombination metrics, it is 
sensible in the 0(0^) calculations we carry out as long 
a-s Rl,r are sufficiently small, e.g. Rl,r < 7r/3 (or for 
any Rl.r using fixed cones). In this section we consider 
different R values to constrain the coefficients of the dif- 
ferent possible NGLs and in the next section calculate 
explicitly the Rl,r dependence of the NGLs for the anti- 
kx algorithm. For C/A and kx, we will calculate for 
Rl ^ Rr = R- Most of the logic and results in this sec- 
tion actually go through for equal i?'s as well. We rely 
on distinct Rl,r only at the very end of this section. 

The dijet events we consider have three kinds of double 
non-local globs, each a ratio of scales in the soft sector: 



g)^a.C./o.(i^.,i?.)ln^^^^^^ 



(£) CF^Af0R{RR,RL)\n^ 
(^yCFCAfMRL,RR)ln'1^ 

\2n/ Kpt 



kr- 



2Atani?fl/2 
fci tani?fl/2 
fltani?L/2 ' 



(17a) 
(17b) 
(17c) 
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The coefRcients /ol, /or, and /lr depend on the jet 
algorithm. The ratios are between the scales (iden- 
tified in [321 133) of soft gluons inside jets fig" = 
fci.i?/ tan(i?i,^ii;/2) and the scale fig = 2A of soft glu- 
ons outside jets cutoff by energy A. Contributions to 
these logs come from three regions of phase space, which 
we label L for the left jet, R for the right jet, and O for 
the out-of-jet region. These regions arc shown in Fig. [2] 
Each NGL gets a contribution from a pair of these re- 
gions, which set the scales in the log. If we consider two 
regions A and B, then the phase space contributions can 
be divided as: 

M{{h})=MA{{h})+MB{{h}) + MAB{{h}), (18) 

where each M imposes a set of measurements on the 
partons in the final state with momenta ki (cf. Eq. ( |B2[ )). 
The first term comes from gluons only in region A, the 
second from gluons only in region B, and the third from 
at least one gluon in both regions. 

RG evolution constrains the relative contributions 
from each of the terms in Eq. (18). NGLs are indepen- 



dent of the renormalization scale but each contribution 
from Eq. ( 18 ) will have /i-dependence. For instance, the 



/or NGL comes from the sum: 
~ (2^-) ^P^^ ■/"oRi^R^ Rl) 



(19) 



2 In 



2 /itani?/j/2 



21n2 4-ln2 
2A 



fj,"^ tan Rr/2 



2Afcf 



The last term is especially notable: it only contains con- 
tributions with two soft gluons in the final state that 
live in separate regions (one in R, one in O), and it 
is the only term of the three that depends on multiple 
scales. These contributions are simpler to compute than 
the other terms, with the added benefit that there are 
no global terms with the same color and log structure. 
These mixed-scale terms alone determine the coefficient 
of the NGLs, as the others are fixed by RG invariance. 
This feature was used in [21] to determine the complete 
set of non-global terms in the hemisphere dijet soft func- 
tion. 

Let us consider the bare contribution to the soft func- 
tion from the last term in Eq. (19 1, following from 
Eq. (B4). To order 1/e^ in the MS scheme, the mixed- 



scale term is 

'-'NG ~ 



(27r)2 r(l - e)2 



X A 



tan 



— 2/oR(i?i?,, Rl) 

2£ 



(20) 



Similarly, for the NGL depending on both and kji, 
the mixed scale term following from Eq. ( B5 ) is 



'-'ng — 



{2n) 



r(i 



2 /lr ( -Rl, i?_R) 



x{kM-'-''tan''^te.n''^ 



(21) 



The 1 /e poles in these contributions are infrared in origin. 
Now, the full soft function is infrared finite. As argued in 
[?rj , this means that the purely "in" and "out" contribu- 
tions in Eq. ( 19 1 contribute compensating IR divergent 



terms that cancel the IR poles in Eqs . ( 20 1 and (21). This 
also cancels the ^-dependent terms in Eqs. (l20| and (pll). 



preserving RG invariance of the factorized cross section. 
After this cancellation, double logs of k^^n/A and kL/kn 
survive in the full soft function, and similar double logs 
of the other scale ratios survive. These are the NGLs. 

The constraints from RG invariance imply relation- 
ships between the NGLs in Eq. (17). It is instructive 



to break up the contributions to the non-global double 
logs in terms of what regions the soft gluons are in. There 
are six such regions, and all of the contributions have a 
coefficient —{as/2i:YCFCA'- 



2 /itan(i?L/2) 



/R(i?i?)ln- 



/o(i?L,i?fl)ln 



2 ^tan(i?fl/2) 



2 JJ^ 
2A ' 



(22) 



-/oL(i?L,i?fl.)ln- 



-foniRR, Rl) In' 



tan(i?L/2) 
2AkL ■ 
/i^ tan(i?fl,/2) 
2Ak^ 



^/LR(i?L,i?fl)ln' 



2 M2tan(i?L/2)tan(i?fl/2) 



Note that the first three coefRcients receive contributions 
from purely real diagrams (with both soft gluons in the 
final state in the same region) and real-virtual diagrams 
(with one soft gluon in the final state and one virtual soft 
gluon). There are several properties of these coefficients: 

• /l = /r , 

• fo{RL,RR) ~ !o{RrtRl) , 

• fm iRL, Rr) — fm iRR: Rl) , 

• fohiRL, Rr) — foR{RR, Rl) , 

• /oL and /or may not be symmetric in their argu- 
ments. 

Finally, the statement that the NGLs are determined 
purely by RG invariance and the mixed scale logs (/lr, 
/ol, and /or) implies that there are relations between 
the coefficients. Expanding the logs in Eq. (17) and using 

Eq. m, 



/l(-Rl) = 2 [/lr (i?L, -Rfl) + foL{RL,RR)] , (23a) 
MRb) = 2 [/lr (i?L, -Rif.) + foRiRR, Rl)] , (23b) 
/o(i?L, Rr) = 2[/oL(i?L, Rr) + /oR(i?fl, Rl)] ■ (23c) 

So far we have not used that Rl r could be different. 



so our proof of Eq. (231 is valid for equal i?'s. Now, we 
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can use different Rl,r to argue that the Rl dependence 
cancels between /lr and /or; similarly, the Rr depen- 
dence cancels between /lr and /ol- If /or is known, 
then up to a constant /lr can be determined. We have 
the additional constraint that as the jet radius shrinks 
to zero, the NGL of the two jet scales will also van- 
ish: /LR(i?L,-Rfl) ^ as i?L ^ or Rr 0. This 
means that only knowing /or completely determines all 
the other coefficients. We will compute /or for the anti- 
kx algorithm and use it to determine /lr. For the C/A 
and kx algorithms we will take Rl = Rr — R- 



promotes a particle to a jet if some di is the smallest. 
This procedure is repeated until all particles are put into 
jets. The common recombination algorithms (kx, C/A, 
and anti-kx) are part of a class parameterized by a real 
number a. In terms of a, the metrics for e^e~ are 

= 2 mm{Ef°',E]"){l - cos %) , 
di = 2Ef"{l^ cos R). (25) 

For both types of algorithms, a veto A on soft jets is 
required for infrared safety of exclusive jet cross sections. 



IV. NON-GLOBAL LOGS FOR SEVERAL 
ALGORITHMS 

In this section we derive results for the non-global part 
5ng of the soft function in Eq. ([7| not predicted by RG. 
For each algorithm we first determine the "in-out" NGLs 
by calculating /or, and then determine the "in- in" NGLs 
by calculating /lr. After calculating the leading NGLs 
for each algorithm, we plot the coefficients of the logs and 
discuss the results. 

As shown in Sec. |III[ the double log terms in ^ng can 
be determined by the calculation of /or, the contribution 
with one gluon in a jet (in this case the right jet) and one 
gluon out of the jets. As is well known, at 0{al) the non- 
global double logs arise from soft gluon emission diagrams 
with the C'fCa color structure, with the amplitude in 
Eq. pSl). 



From the form of Eqs. (B4) and (B5) one finds that 



the coefficients of the leading NGLs are given generically 
by the integral (cf. [10]) 



/o: 



L.ORXR 



= 2 



driidrj; 



cos ( 



TT cosh(77i - 772) 



xeoloR.LR: (24) 



where ryi_2 = lncot(6'i 2/2) are the (pseudo-)rapidities of 
gluons 1,2 with respect to the z axis (the jet 1 axis). 
The angular constraints of the jet algorithm are given in 
©aig, and depend on which coefficient (OL, OR, or LR) 
we are calculating. This integral is a geometric measure 
of the size of the region into which the two soft gluons 
are allowed to go for a given contribution to the NGL. 

We consider two types of jet algorithms, cone and re- 
combination. Cone algorithms fit jets to a geometric 
shape (the cone), and a jet is found when the momen- 
tum in the cone is aligned with its axis. Therefore, the 
phase space constraints for particles in a found jet simply 
requires them to be within an angle R of the cone axis. 
For soft particles in the n jet, for instance, this implies 
dns < R- 

Recombination algorithms build jets by successive 2 — > 
1 mergings of particles. A pairwise metric dij and a single 
particle metric di govern the recombinations. A single 
step in the algorithm finds the smallest of all dij and di, 
then merges the closest pair if some dij is smallest or 



A. Cone or anti-kx algorithms 

To leading order in the SCET power counting, the 
phase space for soft particles that get combined into 
the jets is the same in the cone and anti-kx algorithms 
[551 133], so they will have the same leading NGL (also 
pointed out in [28 ). We work with fixed cones, but the 
results apply to other infrared-safe cone algorithms (e.g. 
27 ) in the configurations we consider where there is no 
split-merge issue. For generic jet configurations, different 
cone algorithms will have different NGLs. 



1. In- Out NGLs 

For these algorithms, the two-particle phase space for 
one soft parton to be inside and one outside of a jet is 



©oT = ^iVR <V1<^) 0{-VL <V2< VR.) 



+ (1^2), 



(26) 



where r/L^n = Incot Rl,r/2. Interchanging the gluons 
simply introduces a factor of 2 into the integral, and so 
we will simply work with the first configuration shown in 
Eq. ( |26[ ) and multiply by 2 to account for this symme- 
try. With these constraints the in-out NGL coefficient 



foRiRR,RL) in Eq. (24) is then given by 



roT{RR,RL)^ dm '^^2 2fa-,.)_T ^ (27) 

Jvr J -til ^ ^ 

The integrand depends only on the difference 771 —772, and 
can be easily integrated to give 

PoT^Rr^ Rl) = ^-2 Li2 Ln^ ^ tan^ ^ . (28) 



we plot Eq. ^ ior Rl 
—T' TT^/S, and at i? -^• 



cone 
OR 

^TiR) ^ 0. 

This result is consistent with that of 



Rr = R. At i? = 0, 
7r/2 (hemisphere jets). 



In F 
/< 

^ in the limit 

i? — > 0, who considered the case of measuring only one 
jet's invariant mass and imposing a jet veto outside the 
two jets. Adding two copies of the NGL they found in 
that case reproduces Eq. ^ (for R ^ 0). Eq. 



now provides the full R dependence of the coefficient of 
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the in-out NGLs in Eqs. (17a I and (17b I for the cone or 



anti-kx algorithms, with the separate Rl,r dependence 
derived for the first tinic."^ 



2. In-In NGL 



We can use the constraints in Sec. IIIII to determine 
the in-in NGL coefficient, fl'^". The sum /or + /lr 
must he Rl independent, and /lr is symmetric in its 
arguments and vanishing as Rl,b. ^ 0. Since /qr''' hi 
Eq. (28 1 happens to be symmetric in Rl.r., it is simple 



to determine 



R 

hniRL, Rr) = 2Li2 (tan^ tan' 



(29) 



This result is simple to check by performing the integral 
in Eq. ( 24 ) with the constraint 



-"LR 



0{m < ??i < oo) 6i(-oo <r]2< til) , (30) 



which confirms the result in Ec^. ( 29 ) 



a constant but noted that no physical explanation was 
apparent. We observe above that the boost invariance 
of the amplitude and the property that all rapidities get 
covered in the sum /£™° + /qr"' removes dependence on 
where precisely that jet boundary 77^' is. Below we will 
find that this coefficient is not a constant for the C/A 
and kx algorithms, where the phase space for recombi- 
nation of two soft gluons in a single jet is distinct from 
anti-kx and changes the value of /r. In particular, the 
phase space included in the sum /lr + /or contains gaps 
not invariant under boosts so that dependence on the 
boundaries remains. 



B. Cambridge- Aachen algorithm 

In computing the leading NGLs for the C/A and kx al- 
gorithms, we will use only one common jet radius, setting 
Rl= Rr = R- In these cases /or = /ol and /r = /l- 



In- Out NGLs 



j-cone/p p 

/lr [Rl,Ri 



-0° 



Rl 

= 2 Li2 ( tan^ tan 



g2(?7i-772) _ I 

2 



(31) 



In Fig. |4] we plot this coefficient ior Rl = Rr — R- 

The sum of the in-out and in-in coefficients is related 
to the coefficient fn, which comes from having gluons 
only in the right jet: 



/•conc/p p \ I fconc/p p \ fco 

/or yn.R,HL) + /lr ytlL.tCR) - -^JR 



\Rl 



(32) 

The fact that /|^°"° is a constant independent of Rr is 
somewhat surprising, but can be shown at the level of 
the integrals in Eqs. (27) and (31). These integrals sum 
simply to give 



df]2 

-00 



16 



g2(r)i->)2) _ 1 



(33) 



The integrand depends only on ?7i — 772 , and we can shift 
the integration variables by 771.2 — ?7i,2 — VRj rendering 
the integral to be a constant: 



a otnc 



dm 



dm 



16 



;2(r,i->)2) _ \ 



2^ 



/o 

In Fis. 5 we compare this constant to the result for Jj^ 



(34) 



from otner algorithms. Ref. 



also noted that /r°"° is 



^ Ref. [19] calculated the same coefficient Eq. ([28} for the NGL of 
the jet veto over the total energy in two-jet events in which the 
jet masses are not measured (for ij^ = Rji = R). 



The phase space constraints from the C/A algorithm 
that figure into /or are more complicated than the cone 
or anti-kx algorithms. To contribute to /or, one of the 
soft gluons must be merged into the right jet while the 
other is not recombined with either jet. This amounts to 
the phase space constraints 



e: 



C/A 
OR 



61(0 < 01 < R) 9{R < 02 < TT ~ R) 

X e{0i < 012) . (35) 



Here 6*12 is the opening angle between gluons 1 and 2, 
while 01 and 02 are the angles between each gluon and the 
right jet axis. The last condition is required to guarantee 
that partons 1 and 2 are not combined together first by 
the algorithm, in which case they would be in the "in- 
in" or "out-out" part of the phase space not contributing 

C /A 

to /or • An analogous constraint only contributes at 
the level of a power correction for the cone or anti-kx 
algorithms. 

The coefficient /qr^ is then determined by the integral 



/o(^(^) 



e 



C/A 



^^'2 / — -OR 
-VR •'0 ^ 

COS 6 



cosh(r7i — 7/2) — cos ( 



(36) 



Here rji^ = lncoti?/2. The phase space constraints in 
®OR^ can be written as 1 minus the region where the 
soft partons 1,2 do get combined by the C/A algorithm, 

f?(ei2-ei) = 1-6(01-012) (37) 
= 1 — (cos0 > niaxje^''^ sinhyyi, —e^^ sinh 772}) . 

The last theta function enforces that partons 1 and 2 are 
closer to each other in angle than the jet axes. The "1" 
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term in Eq. (37) produces the same result as /qr"', so 
that Eq. ( 36 ) can be written as 



C/A, 



(38) 



To evaluate A/q^'^, we can perform the integral over 
using (for r/ > 0) 

rns rh 

del)- 



cosh rj — cos ( 



tanh rj 



tan 



tan(0/2) 



tanh(?7/2) 



(39) 



and evaluate the integrals over 771,2 numerically. The 
result is plotted in Fig. [3] 



2. In-In NGLs 



C/A 

The contribution to /^^ requires one gluon in each 
jet. The two gluons must merge with the jets be- 
fore merging with each other, which requires 612 > 
max(0i , 6*2 ) . Therefore 



e^^"^ = 0{Q < 01 < R) 9{tt ~ R < 02 < tt) 
X 9{max{9i,02) < 612) ■ 



(40) 



Again we can divide up these constraints into a contri- 
bution identical to the cone restrictions and a correction 
factor, 

/lT(^) = flTiR) - A/LT(i?) , (41) 
where the last term is given by 

X ('?2 — lnsinh77i) < — cos^^ (e"''^ sinhTyi) (42) 



-t- 



tanh(77i-772) 



-tan 



coth 



Vi ^ V2\ / 1 ^ 6 sinh 7]i 
2 y Y l+e-''2 sinhT^i 



Considering the phase space constraints, it is straight- 
forward to see that if i? < 7r/3, the in- in NGL coefficients 
must be the same for cone, anti-kx, and C/A (as well as 
kx). This is because for a sufficiently small jet radius, the 
soft gluons cannot merge together since 612 > R- This 
is confirmed by the calculation of /lr^i ^nd is shown in 

Fig.m 

C/A C/A 

With the results for /q^^ and f^^ , we can combine 

C/A 

them to determine . Unlike the anti-kx algorithm, 
this coefficient is R dependent for the C/A algorithm, 
and is plotted in Fig. [5] 



C. kT algorithm 

The phase space constraints from kx algorithm are 
more complex than the C/A algorithm, as there are dif- 
ferent limiting cases to consider where the NGLs can 
be large. To determine whether two soft gluons are re- 
combined by the kx algorithm or not, we consider the 
pairwise distances among the soft particles and the two 
collinear jets, 

di2 - 2iJim{Ef,E^) (1 - COS6I12) , 

d^n = 2Ef{l-cos0.,), d,n = 2Ef {1 + cos 0^), (43) 

for i = 1,2. The single particle metrics are 

d,^2Ef{l-cosR). (44) 

The phase space constraints from the kx algorithm de- 
pend on the relative scaling of the two soft gluons. For 
instance, if gluon 1 is in the R jet and gluon 2 is out 
of both jets, then the limit Ei <C E2 will give rise to a 
different phase space than the opposite limit, E2 <^ Ei. 
Both limits will give rise to a large NGL, and we will con- 
sider both regimes in calculating the coefficients for the 
different logs. The limit Ei ^ E2 will produce the same 
phase space constraints as the C/A algorithm, meaning 
that the NGLs will be the same for kx and C/A in this 
limit. In the limit E2 Ei, we find a new coefficient for 
the NGLs. 



1. In- Out NGLs 

As with the previous algorithms, the coefficient /q^ is 
given by 



dri2 



-IJR 

COS 6 



e 



OR 



(45) 



cosh(77i — 772) — COS0 
The phase space restrictions are given in general by 

©OR = ^(0 <6i<R)0{R<02<T^- R) 

X 0{mm{din,d2) < di2) . (46) 

The first two theta functions require that din < di and 
d2 < {d2n, d2n}, whilc thc last requires that the two soft 
gluons are not merged together by the algorithm. 

The ordering in the last theta function in Eq. ( 46 ) 
depends on the soft gluon kinematics. We consider 
two limits, both of which give a large NGL coefficient: 
kji < Atani?/2 and Atani?/2 < fc^. 

If kii < Atani?/2, then 



i;i(l - cos6'i) <£;2tani?/2. 



(47) 



Unless (?i <C 1 (which is a power correction for the con- 
tribution to the NGL coefficient), then this limit im- 
plies din < d2. This means the relevant comparison is 
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In-Out NGL coefficient 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



R 

FIG. 3: "In-Out" NGLs for three algorithms. The co- 
efficient fo^{R) (equivalently /ql) of the leading NGL 
In^ [fc_R/(A tan 7?/2)] for the cone/anti-kx algorithms (solid), 
the Cambridge-Aachen algorithm (dashed), and the kx algo- 
rithm when kn <^ A (also dashed) and when A <C fc_R (dotted). 
These algorithms recombine soft gluons in successively larger 
regions of phase space, reducing the coefficient of the NGL. 



In-In NGL coefficient 

3.5, , , , 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



R 

FIG. 4: "In-In" NGLs for three algorithms. The co- 
efficient fl'liR) of the leading NGL \n^{kL/kR) for the 
cone/anti-kT algorithms (solid), the Cambridge- Aachen algo- 
rithm (dashed), and the kx algorithm in the limit k^ <^ kn (or 
kn <^ kL) (dotted). The coefficients only differ for R > 7r/3, 
the smallest angle for which the algorithms can begin to re- 
combine soft gluons in separate jets. 



din < di2, which simpUfies to the same constraint as the 
C/A algorithm. 



h < 0i2 



(48) 



In this region the coefficient /qr is the same as for the 
Cambridge/ Aachen algorithm. 

In the opposite limit of the NGL scales, Atani?/2 <C 



E2i&nR/2 < £;i(l - cos6'i) , 



(49) 



Outside of the power suppressed region of phase space, 
^2 < din, and the constraint becomes d2 < di2. This is 
more restrictive than the C/A case, equivalent to 



R < 



'12 ■ 



(50) 



In this region the coefficient of the NGL of A/k^ can be 
written similarly to Eq. pS]), 



(51) 



foliR) = foTiR) - ^foUR) 

where A/q^j^ is given by 

Af^liR) = 4 / dm dm 

COS(j) 



cosh(77i— 772) — COS 
whose (f) integral can be evaluated to give 



9iR-0i2), (52) 



dri2 



log cot[max(7r — 



(53) 



tanh(77i-?72) 



tan 



tan(0o/2) 
tanh 



where 



Vq = cos 



R 



R 



cos — - cosh(77i— 772) — sin — - cosh(77i-|-772) 



(54) 

and 6*1 = 2 cot ^ e^^ . We perform the remaining two 



integrals in Eq. (53) numerically and plot the result in 
Fig.i 



2. In-In NGLs 

When one soft gluon is inside each of the two jets, we 
require 



din 7 d2n < di2 



(55) 



As with /or, we want to consider the limits of a large 
NGL, namely <C k^ and fc^ <C k^. We can con- 
sider these limits simultaneously, taking kn <C k^ for 
definiteness. At leading power this implies Ei <C E2 and 
din *C d2n, since we can neglect the small angle region 
near the jet axes. The constraint therefore becomes 



In < "12 



(56) 



and we have 



= e{O<ei<R)0{7r-R<92 < TT)0{9i<ei2) . (57) 

When k^ ^ kp, the last constraint changes to O2 < O12, 
but this results in an identical coefficient. As before, we 
can divide the calculation as 



(58) 
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Sum of NGL coefficients 
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FIG. 5 

rived in Sec, 



Sum of N GL coefficients for three algorithms. As de- 
the sum /or + /lr is equal to the 
■^GLs from gluons just in the R jet. We 



III 



contribution to 

plot |/]^'^ for the cone/anti-kT (solid), Cambridge- Aachen 
(dashed), and Ict algorithms when kn <^ A (dot-dashed) and 
when A < fcfl (dotted). 1/^°"'' is a constant, 7r^/3. 



where 



riR J-oo JQ TT cosh ( 771 -772)- COS 

X eiOi - 9i2) , (59) 



where 

A/L'S(i?)=4/ dm 



friR »'iniii( — ?7h iln sinh ?7i ) 

X < — cos~^ (e~^^ sinh77i) 



dm 



(60) 



tanh(77i-r/2 



tan 



coth 




1 — e ''2 sinhryi 
l+e^''2 sinhr^i 



We perform the remaining integrals numerically and plot 
the result in Fig. |4] We also sum the in-out and in-in 
coefEcients, /q|j + f^"^ = j/r^i and plot it in Fig.js] 



D. Leading NGLs 

With the leading NGLs in hand, we can compare 
the coefficients between algorithms and learn about the 
structure of the NGLs. These observations confirm and 
extend previous studies of clustering effects on NGLs (e.g. 
pni[2^[5U] ). Our extensive investigation of the algorithm 
and R dependence found above makes clearer this con- 
nection between clustering algorithms and NGLs and al- 
lows us to make more specific statements about the field 
theoretic origin and properties of NGLs. 

We first build a picture of the NGLs that we can use 
to interpret the results for different algorithms. For all 



algorithms, the matrix element contributing to the coef- 
ficients /or, /ol, and /lr is proportional to 



cos A0 
cosh(A77) — cos A(j) ' 



(61) 



where is the azimuthal angle difference and Arj is the 
pseudorapidity difference between the soft gluons. The 
phase space restrictions require the gluons to be in dif- 
ferent regions, so that A77 > 0. The region A77 and A0 
near zero provides a collinear enhancement to the NGL 
coefficients, but there is no collinear singularity that sets 
the value of the coefficient. This implies that the NGLs 
are soft logarithms. The infrared singularities that are 
contained in the distributions 



and A 



-l-2e 



(62) 



come from the soft region of phase space, when £'12 ~^ 0. 
The coefficients of the NGLs receive support over the en- 
tire region of phase space, with the dominant contribu- 
tion near the jet boundary. Since different jet algorithms 
merge nearby soft gluons in different ways, the coefficient 
of the in-out NGLs is very different. We now discuss the 
results for each algorithm. 

We start with the in-out coefficient, fQ^{R) (equiva- 

lently /ql)- Fig- [^j we plot this coefficient as a func- 
tion of R. The cone and anti-kx coefficient is significantly 
larger than either C/A or kx. This is a manifestation of 
the in-out phase space. For the cone and anti-kx algo- 
rithms, the soft gluons on either side of the jet boundary 
will not recombine together. The soft gluons that are 
close to the jet boundary (one in the jet, one out) con- 
tribute more to the NGL from a collinear enhancement 
in the matrix element (see Eq. (61 1). For the C/A and 
kx algorithms, soft gluons near the boundary will recom- 
bine together and not contribute to /or. This removes 
the region with collinear enhancement and subsequently 
reduces the size of the coefficient. Since the kx algorithm 
also weights the pairwise recombination metric by the 
minimum energy of the pair, the region of phase space 
for soft gluon recombination is larger. This further re- 
duces the size of the coefficient, and for large R it even 
changes sign. For each algorithm, it is interesting to see 
that the NGL coefficient is nearly constant over a wide 
range of R; for the anti-kx algorithm, the corrections for 
small R go as i?'*. 

The behavior of the in-in coefficient is very different 
from the in-out coefficient. The in-in coefficient, /lr(^)j 
for each algorithm is plotted in Fig. |4] Each algorithm 
gives the same coefficient for R < 7r/3, and the coeffi- 
cients differ for larger R values. This comes from the 
action of the algorithm: for R < 7r/3, the two gluons are 
separated by more than R, and so they cannot be recom- 
bined. Furthermore, since each jet contains only one soft 
gluon, the action of each algorithm is the same, and so 
/lr(-R) is the same for each algorithm. 

Note that as i? — )■ 0, the phase space shrinks to 
and the coefficient vanishes. This supports the picture 
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that the logs proportional to /lr are soft logarithms, 
since we would expect collinear logarithms to have a 
non- vanishing coefficient in the small R limit. The small 
size of /lr away from the small R limit comes from the 
fact that the soft gluons must be separated by an angle 
greater than tt — 2R. This cuts out the region of phase 
space with a collinear enhancement in the coefficient and 
reduces its size. As R increases, the soft gluons can be 
recombined across the jets, which also reduces the mag- 
nitude of the NGL coefficient. 

Finally, in Fig. [5] we plot the sum of the NGL coef- 
ficients, y^HR) - /o|,(i?) + /lr(^) (or equivalently, 
1 jaig-j^ For the kx algorithm, we have plotted the coef- 
ficient in the two limits A <C Qp and Qp <^ A. For the 
anti-kx algorithm, this coefficient is a constant for all R. 
This comes from the boost invariance properties of the 
matrix element that make it independent of the jet ra- 
dius, see Eq. ([34]). Since /lr is the same for all algorithms 
for R < tt/S and the R dependence of /or is different for 
each algorithm, this implies that /r cannot be a constant 
for C/A or kx. The coefficient /r is the contribution to 
the NGLs from the region of phase space with soft glu- 
ons only in the (right) jet. This coefficient can receive 
contributions from real configurations with two soft glu- 
ons in the final state or real-virtual configurations with 
one soft gluon in the final state and one virtual gluon. 
Therefore the fact that the anti-kx coefficient is constant 
is possibly due to an accidental cancellation. The fact 
that /r coefficient does not vanish as i? — > indicates 
that at least part of the coefficient receives contributions 
from a collinear log. 



V. COMPARISON OF GLOBAL AND 
NON-GLOBAL LOGS TO EVENT2 

We can test for the presence of the NGL of A/mi.2 
in tT(mi, TO2, A) implied by Eq. (17 1 by comparing pre- 
dictions with and without it to the output of EVENT2 
[m 132] ■ To focus on this NGL let us consider measuring 
the two jets' total invariant mass p = {rn\ -\- m\)/Q'^. 
This prevents NGLs of mi/m2 from contributing. But 
first, we must construct the prediction for the global logs 
of p and A/Q in the cross section. 



A. Global Logs 



The resummed cumulant cross section, 



S(p,A) 



: y ' dp' J dmldml S (^p' - 



X / dA'o-(mi, m2. A'), 



(63) 



also splits into pieces predicted by RG evolution and 
those that are not and contain the NGLs, 



I](p, A) = Ei„(p)Eout(A)SNG(p, A) 



(64) 



Predictions for the "global" pieces predicted by RG evo- 
lution, Ein and Sout, are derived in App. [A} 



For the in-jet contribution. 



Si„(p) = 1- 
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(65) 



up to a term constant in p. Here as = as{Q), Lp — Inp, Lji — lntan(i?/2), and ^{R) is the unknown part of the 
0{a^g) non-cusp anomalous dimensions in Eqs. (A9| and (AlO). Meanwhile, the out-of-jct contribution is 



t(A) = l- 
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plus a term constant in A. Here La = ln(2A/(3). In the 
following we only study double logs and L\ , so we do 
not need to know 7(i?). Note that Eout — >■ 1 in the limit 
i? — ?> 7r/2, as it must for hemisphere jets for which there 
is no "out" region. 



The non-global term Sng in Eq. ( 64 ) predicted by the 
results in Sec. IIIII and Sec. IIVI is 



(P,A) 



2 , ^ o2Atanf 



(67) 

with^/^P given for the cone and anti-kT algorithms by 



(36) 



Eq. (27), for the C/A and kx algorithms by Eq. 
Soft gluon clustering, present for the C/A and kx algo- 
rithms, also affects the Abelian terms in the cross 
section [221 [SS]. These clustering effects do not come 
from RG evolution, and can be included as a correction 
factor similar to Sng- 





700 
/ uu 




ouU 




JUU 


dCTEVl 


4UU 


dp 


300 




200 




100 








cone or anti - kj algorithm 




FIG. 6: The CfCa terms in the distribution da^Yi/ dp, using 
the anti-kx algorithm. The coefficient of (jo{ols/2-kYCfCa is 
plotted. Five R values are shown, and A = O.OIQ is used. 



B. Comparison to EVENT2 

For comparison, in Fig. [6] we plot the CfCa terms in 
da^Y2/dp before any terms are subtracted. This gives 
the size of all terms, and after subtracting the global logs 
we can see the size of the non-global terms. 

We can test the predictions for the leading NGLs by 
comparing to the output of EVENT2 [12 [H]. EVENT2 
can numerically calculate an observable vanishing in the 
two-jet limit for e+e~ collisions at 0{a1). For a given R 
and A, we can use EVENT2 to find the p distribution at 
O (a J ) , di7EV2 /dp. By going to the regime p <^ Atani?/2, 
we can numerically enhance the NGLs relative to the non- 
logarithmic non-global terms. 

EVENT2 calculates a binned distribution in p, with 
the cross section in a bin given by 



dp — 1^ = SEV2(Pmax) " SEV2(Pmin) , (68) 
P.„.„ '^P 



with 



Sev2(p) = / dp' 
Jq 



I daEV2 
dp' 



(69) 



Ref. 1301 expresses this NGL as a function of Qp/(2Ai?,^), work- 
ing in the small R limit. The extra factor of R relative to 
Eq. l |67| is due to the imposition in 1301 of a hard cutoff Q/2 
on the soft gluon energy in full QCD. Since the angle is cut 
off by R, this corresponds to cutting off the virtuality of soft 
gluons by ~ Q^tan^(R/2). This choice of cutoff thus cor- 
responds to evaluating the SCET soft function Eq. l |20| at the 
scale /X = Qtan(i?/2). We constructed the global logs in Eq. |64j 
by running the jet and soft functions to /i = Q , which leads nat- 
urally to writing the NGL in the form Eq. ((67|. Thus the NGLs 
quoted in Eq. l |67[ l and 1301 are compatible, with the two conven- 
tions differing in which terms get grouped into the soft NGLs. 
Predictions for physical cross sections remain equivalent. 



In comparing to EVENT2, we will subtract the global 
and leading NGLs: 



Aaip) 



da 



EV2 



dp 



f global 

V dp 



dau 



dp 



(70) 



If only single logs remain, then the resulting binned dis- 
tribution will be constant with value of the single log 
coefficient. The non-log terms in the distribution will be 
a small correction to flatness in the limit p <C Atani?/2. 

For each algorithm, we calculate Aa{p) for A = O.OIQ 
and five values of i?: i? = {0.1,0.6,1.0,1.16,1.3}. R = 
1.16 corresponds to the jet radius used in the EVENT2 
studies in '40'. For each algorithm and R value, we plot 
in Fig.[7]i, Fig. [8^ and Fig.[9 1 the CfCa terms in Ao-giobai 
with only the global logs in Eq. ( 70 ) subtracted out, plot- 
ted in units of cro{as/2n)^CFCA- 



The resulting distribu- 
tions grow linearly for small p, indicating the presence of 
a remaining double log. In Fig. [7)d, Fig. [8]d and Fig. |9]d, 
we plot the CfCa terms in Aa{p) after also subtract- 
ing out the double NGL in Eq. ( [70| . The distribution is 
convincingly flat in the small p regime for all three algo- 
rithms. This means that the double logs of p have been 
successfully removed, confirming the calculations of the 
NGLs performed in Sec. 
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In Fig. 171 we plot the CfCa terms in both A(Tgiobai(/o) 
and Aa{p), where Ao-giobai only has the global logs re- 
moved (and not the leading NGLs). We can clearly see 
the double log dependence in AcTgiobai, and that it is 
entirely removed within statistical uncertainties in Act. 
Note that as R decreases, if we want the non-global dou- 
ble log to be numerically large then the range of p must 
move to smaller values. 

In Figs. Is] and |9j we plot the CfCa terms in Aa{p) 
for the C/A^and kx algorithms. As with anti-kx, we can 
see that the double log dependence has been removed 
within statistical uncertainties. Additionally, we confirm 
the prediction that the leading NGL is the same for the 
C/A and kx algorithms in the limit p <C A. 
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(a) (b) 

FIG. 7: The difference between EVENT2 and (a) the global logs and (b) the global and leading NGLs for the p distribution, 
using the anti-kx algorithm. The coefficient of ctq{ois/2-kY'CfCa in the difference is plotted. Five R values are shown, and 
A = O.OIQ is used. Each difference becomes flat for small p, indicating that only single logs in the distribution remain. 
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FIG. 8; The difference between EVENT2 and (a) the global logs and (b) the global and leading NGLs for the p distribution, 
using the C/A algorithm. The coefficient of aQ{as/2'K)^CFCA in the difference is plotted. Five R values are shown, and 
A = O.OIQ is used. Each difference becomes flat for small p, indicating that only single logs in the distribution remain. 
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FIG. 9: The difference between EVENT2 and (a) the global logs and (b) the global and leading NGLs for the p distribution, 
using the kx algorithm. The coefficient of ao{as/2Tv)'^CFCA in the difference is plotted. Five R values are shown, and A — O.OIQ 
is used. Each difference becomes ffat for small p, indicating that only single logs in the distribution remain. 
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C. Remarks on Soft Factorization 

KSZ proposed that the soft function S'(fci,fcfl, A) m 
Eq. ([6]), at least in the regime A < kL^R <^ QR <C Q, 
factorizes to all orders in as into an in-jet and out-of-jet 
piece, 



S{kL,kR, A) = Sin{kL,kR)Sout{^) , 



(71) 



without any additional factor like S'ng in our Eq. ^ con- 
taining logs of A/kL,R. KSZ recognized that Sin{kL,kR) 
can contain logs of ki^/kR and does not factorize naively, 
but went on to claim that the out-of-jet piece factors off 
completely. Their proof relied on the absence of any logs 
of A/k^^R in S{kL,kR, A). We have shown above that 
such NGLs do in fact arise at 0{a'^) for any value of the 



ratio A/kL,R and that therefore Eq. (71) does not hold. 
Instead 5* takes the form we derived in Eq. ([t]). 

KSZ provided evidence for their claim by comparing 
their prediction, assuming NGLs of A/k^^R are absent, to 
the ©(a^) QCD prediction from EVENT2, using the C/A 
algorithm with R — 1.16 and A = O.OIQ. From our cal- 
culation and EVENT2 analysis, we observe that the C/A 
NGL coefficient for this choice of R and A is too small 
to be resolved in the corresponding plot of HHj (Fig. 1 in 
vl or Fig. 7 in v2), which plots the distribution down to 



p ~ 10 Our calculation of /qr^ in Eq. ( |36[ ) and com- 
parison to EVENT2 shown in Fig. [8] demonstrate that 
the NGL 2a2CFCyi/o/'^(i?)ln2(2Atanf /Qp) is clearly 
present (which, as one realizes from our calculations 
above, is also implicit in the results of [30]). Furthermore, 
based on the results for NGLs in hemisphere soft func- 
tions in [31] we expect single NGLs of A/k^^R and 
non-global non-logarithmic functions of A/k^^R to arise 
in S{kL, kR, A) and also violate the KSZ ansatz Eq. (71 1. 



Our calculation makes clear that NGLs arise from soft 
gluons anywhere in the separated regions of phase space 
into which they are allowed to go, and not only from those 
splitting right along the boundary. Although the latter 
give the largest numerical contribution, as seen in the 
behavior of Eq. (24) as 771 — 772 — > 0, the enhancement 



is not parametrically large, and gluons with significant 
angular separation still give non-negligible contributions 
to NGL. Thus observables sensitive to soft gluons in any 
separated regions of phase space probed with different 
soft scales are prone to NGLs. From the EFT point of 
view, sensitivity of soft modes to any disparate soft scales 
generates non-global dependence on the ratio of those 
scales in the soft function. 



Going beyond the explicit calculations of the leading 
NGLs of jet masses and vetoes, we uncovered strong rela- 
tions among the coefficients of the in-in and in-out NGLs 
and the contributions to NGLs from diagrams with glu- 
ons in the same region of phase space (both gluons in the 
same jet or both outside the jets). We did this by ex- 
tending our insight in [H] that in EFT, NGLs are built 
out of contributions from regions of phase space where 
two gluons enter separated regions or the same region, 
each contributing a log of the factorization scale /i over 
a single scale — k^^R if both gluons enter the same jet, 
A if both lie outside the jets, and intermediate scales 
^/kLjiA or y/khkR if the gluons enter separated regions 
(cf. Eq. (22)). IR safety of the soft function and RG 
invariance require the coefficients of these logs all to be 
related so that they add up to the /i-independent NGLs 
of kL/kR and A/k^^R- By considering the possible de- 
pendence of these coefficients on separate jet radii Rl,r, 
we derived new relations among the coefficients of all the 
different contributions to NGLs. 



VI. CONCLUSIONS 

We have derived the leading NGLs of ratios of jet 
masses rni.2 and an energy veto A on additional jets in di- 
jet cross sections a{m\,m\^ A) using several different jet 
algorithms. We confirm earlier qualitative findings about 
the effects of jet sizes R and clustering in recombination 
algorithms on the size of NGLs involving a jet veto pa- 
rameter, and in addition provide the full algorithmic and 
R dependence for NGLs involving both jet masses and 
vetoes for the first time. Ours is also the first explicit 
calculation in the framework of effective field theory of 
the leading 0{al) NGLs in a soft function appearing in 
the factorization theorem for an exclusive jet cross sec- 
tion involving a jet veto. 

We confirm the insight of [20] that the action of clus- 
tering by recombination algorithms reduces the size of 
NGLs in general. The larger the phase space in which 
soft gluons can be recombined, the smaller the NGLs. 
This is also borne out by the R dependence of the in-out 
NGLs and in-in NGLs. 



The above lessons about the properties of and methods 
to calculate NGLs are directly applicable to exclusive jet 
cross sections in hadron collisions such as at the LHC. 
Gaining control of and resumming NGLs or using meth- 
ods that minimize their impact on jet cross sections will 
be essential to achieving precise theoretical predictions. 

In the quest to resum NGLs using an EFT framework, 
a physical picture is important to recognize non-global 
observables and understand the implication for factor- 
ization theorems. While some elements of the physical 
picture drawn above have previously been noticed, they 
are often not appreciated sufficiently so as to make ob- 
vious the structure of NGLs and when they will appear 
in cross sections. We have studied in detail the NGLs 
that appear in measurements of jet masses in exclusive 
jet cross sections with a jet veto. We hope that this 
work has made clear the source of NGLs, providing both 
better intuition for non-global observables and a quanti- 
tative EFT-based approach to study them. 
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Appendix A: Global Logs from RGE 
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where the constant ji? in Eqs. (A3b| and (A4b) is 2 for 
F = J and 1 otherwise. 

The anomalous dimensions Tf,^f appearing in 
Eqs. (A3) and (A4) are given by 



n=0 n=0 

.(A5) 

where the F f pieces for each function F are proportional 
to the cusp anomalous dimension. 



The parts of the cross section Eq. (|63 
RG evolution are given by (cf. |33l |3l 
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where in Eqs. (Al) and (A2), il 



2uJj{p,Pj) + 

2w™(/i, /is), and each factor Kp = Kp{^, fip),ujp = 
u]f{^,^p) is given by the cusp and non-cusp parts of 
the anomalous dimension of F by 



Kp{f,,f^o)^ r ^ (2rp [a, (a*')] In ^ + 7f [a. {p')]) , 

(A3a) 
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with 
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and the non-cusp anomalous dimensions are given for the 
inclusive jet function by 
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(A8a) 
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for the in-cone soft function by 

1 / 808 11^2 

7s = f^pCA I + — ^ + 28C3 

and for the out-of-cone soft function by 
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(A8b) 

(A9a) 
(A9b) 



, 7o 



2CFrcusp In tan" -w+liR)- (AlO) 



The anomalous dimensions Eqs. (A9) and (AlO) have 
not yet been calculated explicitly beyond 0{as) for jets 
with non-hemisphere radii R, but can be deduced from 
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the known hard and inclusive jet anomalous dimensions 
to 0{ag) and consistency of RG evolution. The relation 
Eq. ( AlO ) between 7out and the universal cusp anomalous 
dimension was proposed in [331 El] based on consistency 
of the RG evolution in the factorization theorem Eq. ^ . 
The term 7(i?) is not yet known, but it must vanish at 
R = 7t/2 and should be non-singular as i? — ?► 0. It must 
also cancel in the sum + 7out • 

The hard anomalous dimensions are constrained by the 
requirement of consistency of RG running to be Th = 



Appendix B: Leading Non-Global Contributions to 
Soft Function at ©(a?) 

The "In-Out" or "In-In" soft functions in Sec. Im] and 



Sec. IV are defined by 



C.OL,OR,LR 
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-2r, 



2r'" and 7ff = -27j - 27^ 



out 
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We also use 



the beta function coefficient /3o = (HC*^ — 2np)/3. 

The fixed-order hard, jet, and soft functions in 
Eqs. (|Al|) and ([A2|) are aU given to ©(a^) by 



X (fci, fc2)A^5°^,f ^'^^1 (fci,^, A)C(fci)C(fc2) 
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where j = {I,T,G,'H,Q} sums over different 0{al) 
cut Feynman diagram topologies (independent emis- 
sion, three-gluon vertex, gluon bubble, ghost bubble, 
and quark bubble, respectively) with two real gluons 
in the final state. The I, T diagrams are illustrated 



Qf 



(All) 



Q in Fig. 10 The cut propagators give factors C(fci^2) = 
^^27r(5(fc2J]F(fcJ 2), where k° = {n-k, + n-k,)/2. The func- 
tions 7\/([OL,OR,LR] jj-Qpose the measurements on the two 
final-state gluons, giving the contributions of the possible 
ways that two gluons can go into separate regions. 



For the jet and in -cone soft functions J, Sin in Eq. ( [AT| ), M[^ JjkL^ji, A) = S{kL-QLn-ki)S{kj^)S{A-Qoutk2 
each log in Eq. (All I should be replaced by the dif- 



ferential operator appearing in the arguments given in 
Eq. (All. The scale Qp appearing in the logs for each 
function F = H,J, Sm , 5'out is 
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Thus it is most natural to evaluate each function F(fip) 
in Eqs. (All and (A2| at the canonical scale iip — Qp, 



+ (1^2), 



(B2c) 



but the expressions Eqs. ( Al I and ( A2 1 are invariant un 



der different choices of fip. 

The constants Cp in the 0{as) fixed-order functions F 
are given by 
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where Ol,r are the phase space constraints restricting 
gluon 1 or 2 to be in the left or right jet, and 9out restricts 
gluon 1 or 2 to be outside both jets. Their form depends 
on the algorithm used to find the jets. 

Expressions for all the matrix elements Aj are given 
in Appendix B of Ref. [21]. In general covariant gauge, 
the leading NGL comes only from the sum of I and T 
diagrams proportional to the color factor CpCa- This 
contribution to the total amplitude is 



(Al3d) 
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n - ki n - ki n - k2 n - k2 ' 
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The one-loop in- and out-of-cone soft function constants 
were first derived in IMj. The two-loop constants c|, 



in Eq. (All) are known for the hard [3S] and jet 



functions, but not yet for the in- and out-of-cone soft 
functions. We do not need the constants cj^ in this paper 
since we only study the logarithmic behavior of S(p, A) 
at 0{al). 



in Z? = 4 — 2e dimensions. Since the matrix elements are 



symmetric in 1 o 2, the symmetrized terms in Eq. (B2) 
just give factors of 2. 

The leading contribution to the non-global parts of the 
soft function S^q in Eq. (Bl) is then given, after inte- 
grating over ki , A:^ and kf^and converting to the MS 
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a) 





b) 




c) 




FIG. 10: 0{al) real emission diagrams contributing to leading NGLs. The endpoints of the gluons can be attached to the 
points on the Wilson lines labeled by a 'x' in any order. Figure (a) gives the I diagrams, (b) and (c) give the T diagrams. 
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For the cone or anti-kx algorithms, the thcta functions 
0i?,L,out take the form 
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H 



2A-k. 
^2 



+ 



(B6b) 



where z = \j k^ A;2^/[fcjj(2A — /c^)]. A similar formula 

holds for S^Q. Meanwhile, the leading contribution to 
SJJq takes the form 



b^Q{kL^R,A) - — r(l-e)2 ^^ 

pOD poo 

X / d/cf (fc/j/cf / dk+{kLkty 



where w= \ k^ k2 / {kRki). 



(B5) 



X — ^ d(j)— -QrQl. 

V7rr(5-e)7o l + w2_2ii;cos0 



Then it b ecomes natural to rescale variables in Eqs. (B4) 

,2 Ri 



and {B5) by k^ = kRx/t&ir^, and fcj = 2 Ay or 
= fciy/tan^ We group the factors of -Rlr and 
2A with the e-dependent prefactors in Eqs. (B4) and 



( |B5[ ) and expand the remaining integrand to leading or- 
der in e, as the subleading terms do not contribute to 
the leading NGL. This results in the particular group 



of terms appearing in S'^q''^^ given in Eqs. (20) and 



(21 1, which we have generalized to other algorithms. We 



also changed variables from the remaining light-cone mo- 
menta in Eqs. (B4| and (B5| to the rapidities of gluons 



1 and 2 relative to jet axis 1 to express the coefficients 
/oR,OL,LR in the form given in Eq. (l24|). 
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